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We present a calculation for rj production in the p— 6 Li fusion near threshold including 
the r\ — 7 Be final state interaction (FSI). We consider the e Li and 7 Be nuclei as a — d 
and a — 3 He clusters respectively. The calculations are done for the lowest states of 
7 Be with J = ( g , ~ J resulting from the L = 1 radial wave function. The 77 — 7 Be 

interaction is incorporated through the n— 7 Be T— matrix, constructed from the medium 
OO . modified matrices for the 77 — 3 He and n — a systems. These medium modified matrices 

are obtained by solving few body equations, where the scattering in nuclear medium is 
taken into account. 
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1. Introduction 

The search for exotic 77— nucleus quasi-bound states «D is guided by the strong and 
attractive nature of the 77 — N interaction in the s— wave^. With this motivation, 
measurements for the 6 Li(p , 77 ) 7 Be reaction^ were carried out by the Turin group 
in 1993 at an incident energy of 683 MeV. A theoretical study of this reaction ^ 
was made and it was summarized that the information available was not sufficient 
to conclude the formation of an 77 — Be quasi-bound state. The interest in this 
reaction has been revived by the recent sophisticated study of this reaction at the 
COSY Laboratory at an incident energy of 673 MeV^. Using the experience gained 
by studying the interaction of 77 with lighter nuclei like the deuteron^ and 3 He^ 
by solving few body equations, we extend our study to a heavier nucleus, i.e., the 
7 Be. 
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Fig. 1. The reaction mechanism for p 6 Li — > r/ 7 Be reaction. 



2. The Formalism 

2.1. The production mechanism 

We consider 6 Li and 7 Be nuclei as a — d and a — 3 He clusters. The reaction 
mechanism is, hence, built by assuming a collision of the beam proton with the 
deuteron in e Li producing an r\ and 3 He. Combining with the spectator a, He 
produces 7 Be in a particular state (see Fig. 1). The production T— matrix for a 
relative angular momentum L between the 3 He and a is written as 



where, Fl{Q) — J °° r 2 dr (ipl Be (r) ]l{Q r) "0o Ll ( r ) ) IS ^ ne transition form factor 
for 6 Li — > 7 Be, with Q = |<f — |fc p . Here, ipQ Lt (r) ancl ^z, Se ( r ) are tne radial 
wave functions for 6 Li and 7 -Be respectively. The T— matrix for the elementary 
process, ( | Tpd-^n^He | ) , is written in a two-step model from our earlier work^A 
The calculations have been carried out by neglecting the effect of Fermi motion on 
( I T pd ^ v z He | ). 

The transition matrix for the p e Li — > r\ 7 Be reaction, including the r\ — 7 Be 
FSI is given as 

T = (k v , m 7 | T p 6 Li ^ n r Be \k p ; m p , m 6 ) (2) 



(q\ T p 6 Li ^ v T Be \ k p ) = i 




M ft 

x ( o ' M I Tpd^ n 3 He I — , m. p ; 1 , m 6 > 



(1) 




(?) w' 7 | T p e Li _, n 7 Be \ k p ;m p , m 6 ) 



(2tt) 3 - + ie 



where fc p and k v are the initial and final momenta in the centre of mass system. 
m p , me and 7717 are the spin projections for the proton, 6 Li and 7 Be respectively. 
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2.2. Final state interaction 

The t] — 7 Be interaction is incorporated through a half-off-shell r\ — 7 Be T~ matrix, 

T v 7 Be (k' kz) = I dxx \ip a iHe(xi)\ 2 {Tz(k' , k, aixi , z)+T a (k' , k , a 2 xi , z)} (3) 

where, T 3 (k' , k , a\X\ , z) and T a (k' , k , 02^1 , z) represent 77 — 3 iJe and 77 — 4 i/e 
scattering matrices, respectively. They are written as, 

t 3 (k' , q, aixi , z) 



1 1 1 f°° 

T 3 (k' , k , aixi , z) = t 3 (k' , k , aiX! , z) + —r / q 2 dq- 

27r Jo 



Hi 

2p 



x T Q (fc' , fc , a 2 a;i , z) (4) 
1 f 00 

T a (fc' , k , a 2 xi , z) = t a (k' , k , a 2 £i , z) + — -5 / q 2 dq ■ 



— / q dq -. 



Z — «— 

x T3(fc' , k , ai^i , z) (5) 

where, x\, the initial Jacobi co-ordinate, is related to the position vector, f\ in the 
He — a — r\ centre of mass system of each of the constituents by r\ — a^Xi. 
z = E — |eo|+iO with |eo| being the energy required for the break up of 7 Be 
into 3 He and a. The t 3 (k' , k , a\X\ , z) and t a (k' , k , a 2 x\ , z) matrices have been 
calculated as in Ref. [7]. 

To begin with, we assume that the ry-meson scatters only once from each of 
the 7 Be constituents, i.e., 3 He and a . The required single scattering terms, the 
t 3 (k' , fc, a\X\ , z) and t a (k' , k , a 2 x\ , z) in Eq. (4 ) and Eq. (5) are constructed 
using the factorized impulse approximation (FIA) | 12 | 13 | j n approximation, 
the momentum in the r) 3 He and r\a. centre of masses are obtained by boosting 
their momenta in the rj — 7 Be centre of mass. This includes the Fermi momenta of 
3 He and a in 7 Be. 



2.3. 6 Li and Be wave functions 

We have used two prescriptions for ^§ L% and ip L Be cluster wave functio ns, nam ely 

(1) Cluster model wave function generated using Wood-Saxon potential SMini and 

(2) Green's function Monte Carlo (GFMC) Variational wave function generated 
using Urbana potential 

Since the energy spacing between the first four low-lying levels of the 7 Be is very 
small, it is important to include the contribution from these levels. In the present 
calculation we consider the lowest states of 7 Be, ^ | ' 5 ) doublet split by 0.48 
MeV, resulting from the L — 1 radial wave function coupled to spin \ of the 3 He. 



3. Results and Discussion 

In Fig. 2, we show the total cross sections as a function of the excess energy. In 
this figure we show the plane wave results and those obtained after including the 
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T] — 7 Be FSI. We find that effect of inclusion of the FSI changes the shape and 
the magnitude of the total cross section curve drastically. In this plot we also show 
sensitivity of the results to different choice of wave functions. The FSI matrix has 
been calculated using an r/N t— matrix which corresponds to a„N = 0.88 + i0.41 
fm. 




Fig. 2. Total cross section as a function of Q for a v n = 0.88 + i0.41 fm. The results are shown 
for two different prescriptions of cluster wave functions using the FIA approach. 




Q (fm 1 ) 



Fig. 3. Form factor, Fl(Q) as a function of Q for L = 1 state for two different prescriptions of 
cluster wave functions. 



In Fig. 3 we show the form factor, Fl(Q), calculated using two different models. 
The region between the dotted lines corresponds to the range of Q, at which Fl(Q) 
in Eq. (1) is getting calculated. This range of Q is shown in Fig. 2 by hashed 
region. It can be seen that two form factors differ in this region which gives rise 
to the difference in the total cross sections calculated using the two different wave 
functions. 

In Fig. 4, we show the angular distributions for three different excess energies 
such that it spans region from close to threshold up to 20 MeV above threshold. 
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Fig. 4. The calculated angular distribution for cluster model wave function using the FIA approach. 



We show results without FSI and with FSI for cluster model wave functions. We 
find that even after including FSI, the angular distribution is isotropic in nature. 
At excess energy ~ 19 MeV, the differential cross section is ~ 1 nb/sr for low-lying 
states of 7 Be resulting from the L = 1 radial wave function, in comparison with 
the experimental number^, i.e., 4.6 ± 3.8 nb/sr. 

A calculation of contributions from the excited states of the 7 Be nucleus to the 
total cross section and that of more detailed FSI is in progress and shall soon be 
reported. 
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